Magnetic Insulator-Induced Proximity Effects in Graphene: 
Spin Filtering and Exchange Splitting Gaps 
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We report on first-principles calculations of spin-dependent properties in graphene induced by 
its interaction with a nearby magnetic insulator (Europium oxide, EuO). The magnetic proximity 
effect results in spin polarization of graphene ty orbitals by up to 24 %, together with large exchange 
splitting bandgap of about 36 meV. The position of the Dirac cone is further shown to depend 
strongly on the graphene-EuO interlayer. These findings point towards the possible engineering of 
spin gating by proximity effect at relatively high temperature, which stands as a hallmark for future 
all-spin information processing technologies. 

PACS numbers: 68.65.Pq, 75.70.Ak, 75.70.Cn, 72.25.-b 



Heat dissipation has become the bottleneck for further 
downsizing CMOS devices, and one of the proposed al- 
ternative is to switch spin without producing charge cur- 
rents. In practice, the fabrication of components able to 
simultaneously inject, manipulate and read out currents 
based on electron spin stands as an overwhelming ma- 
terial and technological challenge. The combination of 
semiconductors with magnetic materials remains to date 
unsuccessful owing to material structural and chemical 
mismatches 

Two-dimensional graphene has demonstrated out- 
standing physical properties such as exceptional electri- 
cal, thermal and mechanical properties p, 0|, but also 
very long spin diffusion lengths up to room tempera- 
ture [7Hl3j. This offers an unprecedented platform for 
the advent of lateral spintronics in which a complete 
integration of spin injection, manipulation and detec- 
tion could lead to ultra-fast electronic circuits compat- 
ible with more-than-Moore CMOS and non- volatile low 
energy MRAM memories [14]. However, a fundamental 
challenge lies in the development of external ways to con- 
trol (gate) the propagation of spin- (polarized) currents 
at room temperature, in view of designing spin logics de- 
vices [m, mil. 

Spin in graphene can be influenced by the presence of 
local magnetic ordering intentionally generated by mate- 
rial design or defects. For instance edge magnetism has 
been shown to develop in graphene nanoribbons (a few 
nanometers wide) for certain edge geometries 03, [Hj, 
or the hole structure of graphene nanomesh [19J was 
also theoretically proposed to offer robust and room 
temperature magnetic states able to affect spin trans- 
port [20] . A lot of interest is also currently devoted to the 
tunability of spin-polarized currents and magnetoresis- 
tance signals by intentional defects, or depositing atoms 
or molecules (such as hydrogen [2l|, [22| or 3d and 5d 



metal atoms [23l-l27j or large molecules [28]). Finally, the 
growth of graphene on magnetic metallic substrates was 
also proposed as a route for tailoring graphene spin prop- 
erties [29l433j . However magnetic conducting substrates, 
which naturally short circuit the graphene layer, restrict 
fundamentally the design of novel types of spin switches. 

In this Letter, we report tunable magnetic proxim- 
ity effects induced on graphene by a nearby magnetic 
insulator. We focus on Europium oxide (EuO) which 
has been recently successfully grown experimentally on 
graphene [HJ. By using ab-initio simulations, with 
both VASP and SIESTA codes, the structure and spin- 
dependent electronic properties of graphene/EuO junc- 
tions are computed. Our findings show that the magnetic 
interaction induces a large spin polarization of graphene 
7r-orbitals, leading about 24% together with a large ex- 
change splitting bandgap of 36 meV. 



The Vienna ab initio simulation package (VASP) [35 



\3l\ is used for structure optimization, where the electron- 
core interactions are described by the projector aug- 
mented wave method for the pseudopotentials [HI, and 
the exchange correlation energy is calculated within the 
generalized gradient approximation of the Perdew-Burke- 
Ernzerhof (PBE) form [H, Eq(. The cutoff energies for 
the plane wave basis set used to expand the Kohn-Sham 
orbitals are 520 eV for all calculations. A 4 x 4 x 1 k- 
point mesh within Monkhorst-Pack scheme is used for the 
Brillouin zone integration. Structural relaxations and to- 
tal energy calculations are performed ensuring that the 
Hellmann-Feynman forces acting on ions are less than 
10" 3 eV/A. Since Eu is a heavy element with atomic 
number of 63, and its outer shell (4f 7 6s 2 ) contains 4/ 
electrons, GGA approach fails to describe the strongly 
correlated localized 4/ electrons of EuO and predicts a 
metallic ground state of EuO, whereas a clear band gap 
is observed in experiment [4l|, |42|. Thus, we introduce 
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FIG. 1: (a) side view and (b) top view of the calculated 
crystalline structures for graphene on top of a six bi-layer 
EuO film, the bottom of EuO is terminated with hydrogen 
atoms, (c) relative energy (to the optimized structure) of 
Graphene/EuO as a function of shifting distance (Az) be- 
tween graphene and substrate. 



a Hubbard-U parameter to describe the strong intra- 
atomic interaction in a screened Hartree-Fock like man- 
ner. For the parameter choice, we fix on-site Coulomb 
repulsion and exchange interaction on Eu 4/ orbital as 
8.3 eV and 0.77 eV, respectively. For oxygen 2p orbitals, 
the on-site Coulomb and exchange parameters are 4.6 eV 
and 1.2 eV, respectively (4l| . 

Using GGA+U method, the structure is first optimized 
from energy considerations and the obtained value of 
5.188 A lattice constant is found to be very close to ex- 
perimental data (5.141 A) with an error of only 0.9% 
and close to the LDA+U (5.158 A). With the optimized 
lattice, we calculated the density of states for EuO with 
ferromagnetic state, where a band gap is observed with 
a value about 1.0 eV. This is consistent with the exper- 
imental optical absorption gaps of 0.9 and 1.2 eV ob- 
served below and above the magnetic transition temper- 
H, El- Here GGA+U method gives better re- 



ature 

suits compared to LDA+U. LDA+U method gives 0.7 
eV band gap with ferromagnetic spin arrangement, and 
1.2 or 1.3 eV for 111 antiferromangetic spin configuration 
(AFMI) or the NiO-type antiferromangetic spin configu- 
ration (AFMII), respectively. 

We next consider the lattice mismatch between 
graphene and EuO. If we use experimental values, 
graphene lattice constant is 2.46 A, and EuO one is 5.141 
A. On EuO (111) substrate, a 2 x 2 surface unit cell is 
about 7.2704 A, which can fit with a 3 x 3 unit cell 
of graphene with lattice mismatch about 1.46%. If the 



GGA+U optimized lattice constant is used, 5.188 A, the 
mismatch is much smaller (less than 1%). In our calcu- 
lations, we used the theoretical lattice constant. 

With such a reasonable lattice matching, we first eval- 
uated the stability of graphene on EuO surface before 
studying spintronic properties. Two structures with 
graphene on oxygen terminated surface and Eu termi- 
nated surfaces of EuO are considered. For these two 
configurations, we have the same amount of atoms, 18 
carbon atoms, 24 oxygen and 24 europium atoms. The 
calculated total energies are found to be -544.50389 eV 
and -545.16824 eV for graphene on O-terminated and 
Eu-terminated EuO surfaces, respectively. One can see 
that with Eu-terminated surface, the system is more sta- 
ble with an energy gain of 0.67 eV. Thus, we use the 
lowest energy configuration of 12 layers of EuO as a sub- 
strate. To avoid the bottom surface effects on graphene, 
the bottom oxygen atoms are terminated by hydrogen to 
simulate graphene on a semi-infinite EuO surface. For all 
calculations, the vacuum length is chosen larger than 20 
A. The optimized distance between EuO substrate and 
graphene is shown in Figure [1] with a vertical distance 
around 2.57 A (nearest C-Eu distance of 2.77 A). 

Using SIESTA package [45[ and the optimized struc- 
ture of graphene on EuO shown in Figure [l(a,b), we 
calculate the local density of states for this system 
[Figure [2(c-h)] with LDA+U for the exchange correla- 
tion functional. The self-consistent calculations are per- 
formed with an energy cutoff of 600 Ry and with a 4 x 
4x1 k-point grid. A linear combination of numerical 
atomic orbitals with double-^ polarization (DZP) basis 
set is used. Due to the existence of EuO substrate, the 
two sublattices of free standing graphene break into six 
folders as shown in Figure [2(a) with different colors. In 
this structure, the calculated magnetic moment of surface 
Eu atoms is found a little bit enhanced, about 7.1 
compared to the bulk value of 6.9 Additionally, the 
sublayer oxygen atoms are found to be spin polarized as 
well, with magnetic moment of about -0.03 \±b- The in- 
teraction with the magnetic substrate remarkably affects 
the magnetic properties of graphene. As shown in Fig- 
ure [2(b), the average spin polarization in the graphene 
layer is found to be about 24%. Here, spin polarization 
is defined as a difference between minority and majority 
states normalized by the total density of states at Fermi 
level, i.e. p = w l~ w ! . 

Since these 18 carbon atoms are broken into 6 sym- 
metry groups, their contribution to the total spin polar- 
ization are different. The spin polarization of magenta 
atoms can reach up to 30.3% [Figure [2(d)], while for the 
yellow and blue ones, spin polarization is about 18% [Fig- 
ure [2(e,g)]. As shown in Figure [2(b), the spin polariza- 
tion is induced on p z orbital, namely tt bond. 

We now scrutinize the proximity effect on the graphene 
band structures. In free standing graphene, the honey- 
comb structure can be seen as a triangular lattice with 
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FIG. 2: (a) The six lattices of graphene on EuO represented with different colors, (b) total density of states of p z orbital of 
graphene, (c)-(h) local density of states on each inequivalent atom of the supercell (a). The spin polarizations in (b) to (h) are 
calculated by comparing the density of states between minority and majority states normalized by the total density of states 
at Fermi level, i.e. p 
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FIG. 3: Band structure of graphene on EuO. Black and red 
lines represent spin up and spin down bands, respectively. 
Inset: zoom around the Dirac cone, the symbols correspond 
to DFT data while the lines correspond to the fit of Eq. J2]). 
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FIG. 4: Band structures for graphene on EuO with graphene 
shifted (a) inward (compared to optimized structure) 0.5 A, 
(b) outward 1.0 A, (c) outward 2.0 A and (d) outward 5.0 A, 
respectively. 



a basis of two atoms per unit cell, with 2D lattice vec- 
tors A = 4f (>/3, 1) and B = 4f(>/3, -1), where a is 
the graphene lattice. Of particular importance for the 
physics of graphene are the two K and K' points at 
the inequivalent corners of the graphene Brillouin zone 
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-^), which are called 



Dirac points. The band dispersion close to the K (or K / 
vector, as k = K + q, with |q| <C |K|, has the form, [4 



E±(q)^±v F \q\ + 0[(q/K) 2 ] 



(1) 



where q is the momentum measured relatively to the 
Dirac point and vp is the Fermi velocity, given by vp = 
y/3tao/2, with a value vp — 1 x 10 6 m/s, and t is the 
nearest-neighbor hopping energy of 2.8 eV. 

In the usual case with e(q) = q 2 /(2m), where m is the 
electron mass, the velocity, v = k/m = ^/2E/m, changes 
with energy. In Equation [TJ the Fermi velocity does not 
depend on the energy or momentum [Hj . 

In case of graphene on magnetic insulator, EuO, as 
discussed earlier, the two sublattices have been broken 
into six groups [Figure [2(a)] and there is spin injection 
to graphene from EuO substrate. The linear dispersion 
of graphene band structure is modified [Figure [3] , with 
a band gap opening at Dirac point. More interestingly, 
this degeneracy lifting at the Dirac point is spin depen- 
dent: we have fitted the band structure with a simple 
spin dependent Dirac dispersion relation in presence of a 
spin dependent mass (gap) term 

E a (q) = ±V(hv a q)* + (A a /2)i (2) 

and obtain gap widths with values of = 98 meV and 
=134 meV for minority and majority states, respec- 
tively while the Fermi velocities v± = 1.40 x 10 6 m-s _1 and 
= 1.15 x 10 6 m-5 _1 are also polarized [see inset of Fig- 
ure [3]. The corresponding polarization, around 20% for 
both gaps and velocities, is very significant. There is in 
particular a large energy window (inside the gap region) 
where the graphene would be 100% polarized (half-metal) 
along the majority or minority direction depending on 
the position of the Fermi level. The observed spin split- 
ting is due to the interaction between C-p z and Eu-4/ 
states. Indeed, there is a strong peak of polarized Eu-4/ 
state right below Fermi level. These polarized states get 
hybridized with graphene, hence the induced magnetism. 
This can also be seen from the band structure, where the 
graphene bands start overlapping with the majority Eu- 
4/ bands. 

Since the interaction between graphene and the sub- 
strate is quite weak (as evidenced by the large equilib- 
rium distance of 2.57 A between graphene and EuO lay- 
ers), it can be further easily affected by the external envi- 
ronment. To mimic the situation of internal pressure and 
strain, we calculated the electronic properties for varying 
the interlayer spacing (Az) as described in Figure [1(a) 
and (c). It is found that by interlayer displacement of less 
than 1 A, the total energy of the bilayer system changes 
only by 0.085 and 0.156 eV/cell for displacement in- and 
outward of 0.5 A, respectively. This weak modification 
(for supercell contains 18 carbon atoms) corresponds to 
fluctuations in energy per carbon atom in the order of 
few meV. 

Even though the energetics stability with interlayer 
spacing Az is weakly affected, the impact on the band 
dispersion of graphene is markedly strong, as seen in Fig- 
ure |U When compressing the bilayer by 0.5 A, more elec- 
trons (and spins) are transfered to the graphene layer due 



to enhanced overlap between C p z and Eu 4/ orbitals. 
Accordingly, the Dirac point is moved deeper inside the 
valence bands compared to the equilibrium situation (cf. 
Figs. [3] and a)). In contrast, for larger layers separation, 
the Dirac cone is clearly seen to be shifted out from the 
valence band of EuO, approaching the Fermi level of the 
system [Fig.|4^b) and (c)]. Simultaneously, with the shift- 
ing of the Dirac point out of EuO valence band, the gap 
between spin-up and spin-down bands is continuously re- 
duced. Finally, for Az = 5 A, the spin up and spin down 
branches become almost degenerated and the Dirac point 
crosses the Fermi level, i.e. approaching typical band 
structure characteristics of isolated graphene [Fig. 5(d)]. 

In conclusion, we have reported first-principles simula- 
tions showing that the proximity of a magnetic insulator 
will induce a strong spin polarization of graphene tt or- 
bitals. The Europium oxide substrate was found to break 
the bipartite lattice of graphene into six inequivalent sub- 
lattices, causing variable spin polarizations on the new 
graphene sublattices with an average spin polarization 
about 24%. Simultaneously, a band gap develops at the 
Dirac point with a large exchange splitting of more than 
30 meV, larger than anticipated [48|. These theoretical 
findings deserve further experimental demonstration of 
spin filtering effect and spin-dependent gap in graphene 
based structures. 

Note, after the submission of this paper, we just no- 
ticed that a successful experimental fabrication of EuO 
on Graphene (49| . 
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